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tice with spacing b to regulate the gauge fields in standard fashion, 
while computing the chiral fermion determinant on a finer lattice with 
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1 Introduction 



Particle theories which go beyond the Standard Model are likely to depend 
on the non-perturbative dynamics of chiral gauge theories. It is well appre- 
ciated that an important step for theoretical and computer investigations of 
these theories is their formulation on the lattice. The major obstacle to ac- 
complishing this goal is the seeming difficulty of avoiding 'fermion doubling', 
while at the same time maintaining chiral gauge invariance . In this paper 
we describe a non-perturbative lattice construction, which we show leads in 
the continuum limit to Green functions satisfying the Ward identities of a 
given chiral gauge theory, precisely when gauge anomalies cancel among the 
fermion representations. 

Several of the tools we employ exist in some form in the earlier literature 
on the subject. The most important of these is the idea of coupling contin- 
uum fermions to interpolated lattice gauge fields as a step towards main- 
taining chiral symmetry. Reference Q contains early work on this topic. 
More recently it has been used in the context of maintaining ungauged chi- 
ral symmetry by 't Hooft |J. Also, the authors of § and very recently ||] 
and H have used a gauge field interpolation to the continuum in defining a 
fermion determinant for chiral gauge theories. 

Our chiral gauge theory construction makes use of two lattices: the gauge 
fields live on the links of a regular lattice with spacing 6, while the fermions 
live on a much finer lattice subdividing the 6-lattice, with spacing / < i. 
The gauge field action is the standard one. For the fermion-gauge boson 
interaction we gauge-covariantly interpolate the Wattice gauge fields to /- 
lattice gauge fields and couple these to the fermions. If the interpolation is 
'smooth' enough this will result in a regulated theory in which the gauge 
boson momenta are cutoff above 1/6, while the fermion momenta are at most 
of order 1//. The utility of having these two cutoffs is easy to understand 
on very general grounds and we outline these next. 

The Nielsen-Ninomiya theorem |?J strongly suggests that in order to re- 
move unwanted doublers from the spectrum of light fermion modes on the 
lattice, the chiral gauge symmetry must be broken explicitly. Resigning our- 
selves to this conclusion, there is a simple lattice fermion regulator which 
eliminates the doubler-modes with a Wilson term and results in a gauge 
non-invariant fermion determinant. However, this breakage of gauge sym- 
metry can be cured by renormalization. The approach of restoring chiral 
symmetries through the addition of counterterms is well-known to work in 
lattice QCD [||. The case of chiral gauge theories has been considered by 



1 



the Rome group M |l(J, and also by the authors of O]. They propose to 
recover a gauge invariant theory by a non-perturbative tuning of countert- 
erms. The major advantage of our approach is that the use of two cutoffs 
allows us to recover a symmetric continuum limit after an analytic one-loop 
subtraction. 

Due to the Wilson term, the fermionic one-loop effective action is not 
gauge invariant. Under a gauge transformation on the lattice, it trans- 
forms as, 

^r = L» + o(/ 2 ), (i) 

where L a {x) is a sum of local, dimension < 4 operators made from the gauge 
fields, while the 0(f 2 ) terms may be non-local in the gauge fields. This form 
is not surprising. Any sensible regulator with short-distance cutoff / should 
give gauge invariant results in the / — > limit, except in the presence 
of ultraviolet divergences, which are sensitive to the details of the cutoff 
procedure. The local nature of these ultraviolet divergences underlies the 
locality of L a . Any other violation of gauge symmetry must vanish as / — > 0. 

The fact that the large violations of gauge invariance are local suggests 
that they can be subtracted by defining a 'renormalized' fermion effective 
action, 

IV = To — S c ,t., 

5 a x S c .t. = L a (x), (2) 

where S c .t. is local in the gauge fields. The reason that the locality of the 
subtractions is important is that the non-local 'finite' parts represent the 
consequences of one- loop unitarity and must be left intact. Clearly T r is 
now gauge invariant up to order f 2 . This procedure provides a nice under- 
standing of chiral gauge anomalies. The (consistent) anomaly [^] is the only 
term in L a which is not the gauge variation of a local functional. Therefore 
it cannot be subtracted. Gauge invariance up to 0(/ 2 ) can only be restored 
in theories for which the gauge anomaly cancels among the chiral fermion 
representations. Obviously, we will assume this to be the case. 

If we had exactly cancelled all gauge non-invariance we could proceed 
to integrate over gauge fields without further anxiety, since the gauge fields 
are gauge-invariantly regulated on the lattice. But this is not the case, 
since there remains an 0(/ 2 ) non-local violation of gauge invariance. When 
the fermion determinant is inserted in gauge boson loops, new divergences 
threaten to overcome the 0(f 2 ) suppression. This is where having two 
cutoffs makes a crucial difference. The point is that the new gauge boson 
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ultraviolet divergences are at most given by powers of 1/6 <C 1//. Therefore 
the violations of gauge invariance are suppressed by the small parameter 
f 2 /b 2 to all orders! This is the central mechanism in our construction for 
reproducing a gauge invariant theory in the continuum limit. 

A regulator with two cutoffs, / <C b, may seem strange at first, but it 
is simple to understand if we imagine integrating down to a theory with a 
single cutoff, b. One could in principle integrate out fermion loop momenta 
between the scales 1// and 1/b, to yield a theory cut off at b, but which 
contains a series of higher dimensional local interactions resulting from this 
integration. The new theory is equivalent to the two-cutoff theory for am- 
plitudes with external momenta below 1/b. In particular, though gauge 
invariance is not manifest, the higher-dimensional interactions will compen- 
sate the gauge non-invariance of the cutoff procedure to yield amplitudes 
related by the gauge Ward identities to order f 2 /b 2 . 

The remainder of this paper is organized as follows. In section 2, we 
begin with the formal path integral for a chiral gauge theory and explain 
how the various parts are to be formulated on two lattices, focussing mostly 
on the fermions. In particular we explain how the loss of gauge invariance 
in the regulated fermion determinant can be compensated by subtraction of 
local counterterms. The obstruction provided by gauge anomalies to doing 
this is emphasized. We also propose a hybrid of this fermion regulator and 
the standard Wilson regulator for vector-like theories. It separately regu- 
lates the real and imaginary parts of the fermion effective action, and enjoys 
improved gauge invariance properties over the original regulator, which will 
be important in section 4. Section 3 describes our adaptation of 't Hooft's 
gauge field interpolation scheme, stressing the important properties it satis- 
fies. In section 4, we study the scaling of the effective action with /, before 
integration over gauge fields, taking into account the / dependence intro- 
duced by the interpolation. We show that the effective action in the hybrid 
scheme is gauge invariant to 0(f 2 ) and finite in the limit / — > 0. In section 
5, using the results of section 4, it is shown that the gauge Ward identi- 
ties (or BRST identities if gauge-fixing is employed) hold to 0(f 2 /b 2 ) after 
functional integration over gauge fields. The choice of how small to take / jb 
in practice is then made. Section 6 deals with anomalies in global currents 
which occur in chiral theories. The Adler-Bardeen theorem is demonstrated 
to hold at one loop and then proven to hold to all orders. The second step 
is made particularly simple by the use of the two-cutoff regulator. Section 7 
outlines how the anomalies in global currents translate into the phenomenon 
of fermion number violation in our construction. Section 8 discusses some 



3 



new ideas for improving computational efficiency. Section 9 provides our 
conclusions. 

2 Two-Cutoff Lattice Chiral Gauge Theory 

Our target theory is one describing a left-handed fermion with a non-abelian 
gauge charge. The continuum version which we want to regulate with the 
lattice is formally given by the generating functional 

e W[J„,rj,n] _ J DA^ThpVtp exp - / z trF ^+^D^+J^+rji>+fr] 

= J VA^ e-JI^+^A, e r[A] e -ffnG[A]^ (3) 

where 

e r ^ = detiD, 
G[A] = (iD)~\ 

t) = (ft + i A^) L + ft R, (4) 

and R, L = |(1 ± 75). We see that tpL transforms (in general reducibly) 
under the gauge group, while i/jr is taken to be sterile. 

The above manipulations are valid assuming that iD has no zero-modes, 
which would leave G[A] ill-defined. However, the occurrence of zero modes 
in topologically non-trivial gauge field backgrounds is physically important, 
for example in fermion-number violating processes. For most of this paper 
we will restrict our attention to fermion number conserving processes. In 
sections 6 and 7 we show that fermion number violation occurs in our con- 
struction and we describe a procedure for defining general fermion number 
violating amplitudes. 

Our lattice construction of the target theory ([|) will take the form 

^ s 

where = e lbA)1 ^ is the gauge field variable on the link (s, s + jl) of a 

regular lattice with spacing b. The gauge boson kinetic terms are defined in 
the usual way on the 6-lattice, 




= U^s)U v {s + fi)U; l {s + v)U-\s). (6) 
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To define the fermion determinant however, we employ a much finer 
lattice which subdivides the 6-lattice, with spacing / <C 5, where b/ f is an 
integer. The fermion determinant is then given by the /-lattice path integral 
with action 

Sfm = y T,(*( n h»l R + « M L]*(n + A) - *(n + ^[R + u- l L\V{n) 

n,fi 

f 3 

+4r/ 3 ^ f (n)tf (n) - ^y 5^[*(n)*(n + A) + + A)* (»)]}, (7) 

where n denotes the points of the /-lattice (as opposed to s for the points of 
the 6-lattice). We refer to this as the chiral Wilson (xw) lattice regulator to 
distinguish it from the standard Wilson (sw) regulator for vector-like gauge 
theories [[0|. It is also the regulator used by the Rome group for / = b. 
The action describes a left-handed fermion coupled to a gauge field u^, and 
a sterile right-handed fermion. The two chiralities are coupled together by 
a Wilson term proportional to r, which gives the left- and right-handed 
doubler modes a Dirac mass of order r/f but clearly breaks chiral gauge 
symmetry explicitly. We take r = 0(1). The explicit breaking of gauge 
invariance is a nuisance but not fatal. As we will show, it can be made small 
to all orders by the addition of gauge field counterterms determined at one 
(fermion) loop. Note that we have not added a mass for the fermion. It is not 
required because the x w action possesses a 'shift' symmetry, ipR — ► ifiji + £, 



where £ is a constant spinor [14]. On the other hand, this is not a symmetry 
of the standard Wilson term of vector-like theories and, as is well-known, 
mass counterterms are required in that context^]. 

The /-lattice link variables, u^{n), are interpolations of the 6-lattice 
gauge fields, U u (s), which are the real dynamical variables. The one-loop 
fermion effective action (the logarithm of the fermion determinant) will 
therefore be a functional of the 6-lattice gauge fields, 

F[U}=T[u[U]]. (8) 

The interpolation u[U] must be carefully chosen so that a (nearly) gauge 
invariant, lattice rotation and translation invariant, determinant functional 
of Un induces a (nearly) gauge invariant, lattice rotation and translation 
invariant functional of £7„. Furthermore, the non-locality of T[u] correctly 

tin any case, in the two-cutoff theory only a mass of 0(f /b 2 ) would be generated. 
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describes the effects of virtual massless fermions dictated by one-loop uni- 
tarity, and it is essential that this same non- locality be found in the T[U] 
induced by the interpolation. In section 3, we will describe such an inter- 
polation and its properties in detail. For the remainder of this section we 
will assume that we have such an interpolation procedure and concentrate 
on compensating the breakage of gauge invariance (in terms of u^) inherent 
to the x w regulator. 

We will denote the unsubtracted (gauge non-invariant) x w determinant 
by expTofii]. Consider the effect of an infinitesimal gauge transformation, 

Ufi(n) -> cj(n)n M (n)u;~ 1 (n + p), (9) 

with u>(n) = e l6,a ( n ) Ta . We have then, 

<Lr [u] = < 5 u sf° > u , 

8 u Sf = -*4/ 3 £ 0»{*(n) 75 T a ^{n) 

n 

+ ^p(n-/)) + 0(n + £)) LT a *(n) + h.c.}} (10) 

where only the Wilson terms have contributed to the non-invariance of the 
action, and < .. > u denotes the fermionic expectation in the background of 
the gauge field. 

Naively the breakage of gauge invariance is only of order f 2 . This is 
because the Wilson term corresponds to a dimension-five interaction sup- 
pressed by / for the physical modes, and to a mass of 0(1/ f) for the doubler 
modes, whose effects one might expect to be suppressed by their mass. Fur- 
thermore, these interactions are chirality flipping and must therefore come 
in pairs. This naive expectation is however upset by power-like ultraviolet 
divergences in 1// in the one loop computation of ([l0[), which can cancel 
the f 2 suppression. Of course these ultraviolet divergent parts are local 
functionals of the gauge fields and u, the allowed forms of which can be 
determined by power-counting, lattice rotational symmetry and invariance 
under space-time independent gauge transformations. The finite parts of 
<5uToM can be non-local in u M , but they are suppressed by the power of f 2 
just discussed. Thus, 

5 UJ T [u]=J2Hn) + O(f 2 ), (11) 

n 

where L(n) is a local function of u and u^. L(n) can contain local operators 
up to dimension 5. 
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We are now in a position to see the difference between the breakage 
of gauge invariance caused by a non-invariant regulator and that caused 
by having an uncancelled gauge anomaly. For example, one of the terms 
allowed in L(n) is 

L i( n ) = ^^^{"{nW{n)u- 1 {n + ji) - u^n)] + h.c. (12) 

It is invariant under lattice rotations and it vanishes for constant uj. The 
power of 1// is dictated by dimensional analysis since / is the only scale 
in Tq. The coefficient a however, must be extracted from a one-loop lattice 
diagram. Although the breakage of gauge invariance given by this term is 
large, the fact that it is local means we can substract it from Tq if there 
exist a local functional of u^, S c ,t.[u] such that S^S^t. = L±. For the present 
example such a subtraction can be performed, 

T r [u] = T -Sl t [u], 

j ii 

which contains a gauge boson mass counterterm. We are only allowed to 
subtract local functionals because the non-local part of Tq already correctly 
reflects the requirements of one-loop unitarity. If a term in L(n) can be 
cancelled by subtraction of a local functional up to 0(/ 2 ) we will refer to it 
as 'locally integrable'. 

In the absence of gauge anomalies, all the possible terms in L[n) are 
locally integrable. This result for the \ w action was first proved in |J at 
one loop. The full set of gauge field counterterms aref| 

• Mass Counterterms: 

f&i 2>[$») - l){u-\n) - 1)] (14) 
J n 

• Wave Function Renormalization: 

Y^Z w Tr ^(A^O 2 + 5(Tr (£, Vv)' + SZ NC Tr ^(A^^) 2 ] (15) 

^We follow the notation of ref. We however do not include the operators involving 
ghosts since they are not needed in our two-cutoff construction. 
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• Gauge Boson Self Interaction: 



SZ A , A ^TrJ2(^ ~ l) 2 rr^(^ - l) 2 , 

SZ 4 , NC ^J2 Tr ^ - l ? K " I) 2 ]- (16) 

The A M are lattice derivatives. So, if L(n) is locally integrable, a one-loop 
subtraction yields 5 w r r [u] = 0(f 2 ) and, despite the lack of a gauge invariant 
regulator, gauge invariance is recovered in the continuum limit / — > 0, at 
least at one loop. 

The only term that can occur in L(n) which is not locally integrable is 
the consistent anomaly, 

L a (n) = -j^2 Tl °( n ) e a/3jS AalapAjas - ^apa^as] + 0(/ 2 ), (17) 

where u^(n) = e l ^ a ^ n \u{n) = e ie ^ a \ A brief description of the derivation 
of this expression on the lattice is given in section 6. The gauge color 
trace is proportional to the anomaly coefficient, d a \, c = ^Tr[T a {T b ,T c }], 
of the fermion representations. If it is not zero, the anomaly breaks gauge 
invariance and cannot be subtracted away with any local counterterms. We 
will assume that all gauge anomalies cancel. Hence T r [u] is gauge invariant 
up to O(f). 

Having dealt with gauge invariance (at one loop), let us consider the 
physics that the fermion effective action incorporates. We have used the 
freedom one has in regulating fields in the ultraviolet to impose two cutoffs 
on the theory. The physics at distances S> b 3> / which survive the con- 
tinuum limit will be insensitive to whether we use one or two cutoffs. In 
principle we could first integrate out all fermion momenta above 1/6 and 
get a local theory with just one momentum cutoff. It would contain extra 
higher dimension interactions dependent on b and /, some of which would 
break gauge invariance. But this breakage would be cancelled by the non- 
invariance of the cutoff procedure to yield just an 0(f 2 ) breaking of gauge 
invariance in the one-loop amplitudes of the theory. 
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This brings us to an important point. We have just seen that all un- 
suppressed gauge non-invariance at one loop can be subtracted away (when 
gauge anomalies cancel) . What remains of the gauge non-invariance is insen- 
sitive to the / — > limit. But there will be one gauge invariant interaction 
induced by integrating out fermion momenta between 1/f and 1/6, which 
will be sensitive to / — > 0. This is the gauge field kinetic term which must 
be renormalized by integrating out fermion momenta between 1/f and 1/6. 
The result for its coefficient is, 

l/flg = l/flg + Iog(6//)^. (18) 

where ti (F) is the sum of quadratic Casimir invariants for the fermion repre- 
sentations. Recall that 1/<7q is the original gauge coupling of the two-cutoff 
theory, multiplying the gauge kinetic term, S g . The second term on the 
right-hand side is just the one- loop renormalization induced by integrating 
out fermion momenta between the two cutoffs. Therefore gt, is to be identi- 
fied as the single true coupling of the effective one-cutoff theory which must 
be scaled to zero in the continuum limit 6^0. In terms of gt,, the action of 
the effective one-cutoff theory is insensitive to / — > (assuming, as we have 
up till now, that our interpolation procedure does not reintroduce sensitivity 
to 1/f). Having identified the correct physics within the effective one-cutoff 
theory, we can proceed to always use the explicitly constructed two-cutoff 
theory for practical purposes. 

In a straightforward one-cutoff theory with cutoff 1/f, the one- loop 
suppression of gauge invar iance breaking by f 2 would not have bought us 
much. The problem is that the real suppression factor is f 2 x gauge boson 
momenta 2 . Once T r [u] is put into gauge boson loops, gauge boson momenta 
of order 1/f could arise, cancelling the f 2 suppression and leading to large 
violations of gauge invariance. This is not the case in a two-cutoff scheme 
such as ours. In loops, the gauge boson momenta are of order 1/6 at most 
since they live on the 6-lattice links, so that deviations from gauge invari- 
ance will be suppressed by 0(f 2 /b 2 ), so long as the interpolation procedure 
is smooth enough to not reintroduce substantial 0(1/ f) momentum com- 
ponents into the interpolated gauge field u (ie. the 6-lattice softly cuts off 
the momenta above 1/6). We can then choose f/b small enough for any 
given precision and space-time volume we wish to work with. Unfortunately, 
we will show in section 4 that the interpolation procedure of section 3 is not 
smooth enough for the x w scheme outlined here to work! We are therefore 
led in the next subsection to a scheme which is a simple hybrid of our \ w 
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scheme and the standard Wilson (sw) scheme for vector-like gauge theories. 
We will show that our interpolation procedure is smooth enough for this 
improved hybrid scheme to work. It is also quite possible, that a smoother 
interpolation scheme can be found, in which case the x w scheme would be 
feasible. 

Finally, we define the fermion propagator in an external gauge field, G[u], 
needed in order to define Greens functions with external fermion lines, 

If we only consider processes without external sterile fermions, then we must 
have an even number of chirality flips from the Wilson term in the propaga- 
tor. Since it is the only term breaking gauge invariance and we are not yet 
integrating over gauge loops, it is easy to see that the propagator is gauge 
covariant up to order f 2 . 

2.1 Hybrid Regulator 

We describe an alternative lattice construction of the one-loop fermionic ef- 
fective action which separately regulates the real and imaginary parts. It will 
automatically be gauge invariant in the limit / — > 0, without subtractions, 
precisely when gauge anomalies cancel among the fermion representations. 
This type of separation was introduced in ref.[|l5|] and further discussed in 
the lattice context in refs.[||[16]. Since the trick used may be unfamiliar we 



end the section with a brief discussion of the legitimacy of the regulator in 
preserving the unitarity property of the unregulated fermion determinant. 

The basic, well-known observation is that formally the magnitude of the 
fermion determinant, e Rer , is parity-even and is just the square root of a 
vector- like determinant, which we know how to formulate with exact gauge 
invariance on the lattice. The phase of the chiral fermion determinant, IuT, 
is parity-odd but can be regulated by means similar to those of the previous 
section, though of course this leads to some breakage of gauge invariance. 
The advantage here is that this breakage is small. As we saw in the last 
section, it is given by the sum of local dimension < 4 terms in the gauge 
field plus non-local terms which are however suppressed by powers of f 2 . It 
is easy to see from the possible local terms listed in the last section that the 
only parity-odd term is the consistent anomaly. For anomaly-free theories 
the breakage of gauge invariance is thus of order f 2 . 
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To effect the separation into Rer and Imr, we note that formaliy, 



it) = [it) {it))^ 1 (iD) f %b]h. (20) 
Taking the determinant of both sides of the above equation yields 

detfiD] = {det[W] detftiD)*] det[iD]}l. (21) 

det [iD]* 

Now, since (iDy = if)R + i$L, while it) = itf)L + i<ftR, the third and fourth 
determinant factors in the square root together describe a Dirac fermion 
with vector-like gauge couplings, and a sterile Dirac fermion. Thus we have, 

det[id] = { ^f^i det[iP] det[« fl}*. (22) 
det [iD)\* 

The last, gauge field independent determinant only contributes a constant 
to the fermion effective action and can be dropped. The first ratio is a phase 
factor, while det[i TJ)\ is a real, vector-like determinant. Its square-root can 
consistently be taken as positive as long as the theory does not suffer from 
Witten anomalies 17]. The idea of this section is to lattice-regulate these 



determinants separately, in a way which manifestly respects the even and 
odd parities of Rer and Imr respectively. 

We now have at our disposal a simple means of formulating our new 
improved chiral fermion determinant, 

det[it)} hybrid = ( f^H i (det[i 0U)i (23) 
det[iD\^ w 

Here, det[il)] x ^ is just the /-lattice chiral determinant of the last sec- 
tion, while det [i-PJstu is the standard gauge-invariant Wilson determinant 
for vector- like fermions. The only difference with single-lattice vector- like 
theories is that a mass counterterm is still not required, because the mass 
induced by breaking chiral symmetry via the standard Wilson term is at 
most 0(f/b 2 ). 

As explained earlier in this section, the only breakage of gauge invari- 
ance in dek[iD]hybrid must be parity-odd, and therefore suppressed by f 2 in 
anomaly free models. 

Let us view the above decomposition and regularization in diagrammatic 
terms. The loop diagrams representing the chiral fermion effective action 
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have vertices containing L = (1— 7s)/2. Even before regulating the diagrams 
one can separate the real (parity-even) and imaginary (parity-odd) parts of 
the original diagrams. It is easy to check that the diagrams representing (the 
unregulated) detfi^?] 1 / 2 and (det[iD]/ detfi-D]*) 1 / 2 reproduce exactly these 
parts respectively. By regulating these determinants separately we are just 
regulating the parity-even and odd diagrams separately. 

This scheme does not harm the unitarity properties possessed by the 
unregulated diagrams. The point is that any diagram can be written as a 
non-local part which is insensitive to the cutoff and a local part (polynomial 
in external momenta) which may be sensitive to the cutoff. Unitarity con- 
strains only the non-local parts, the local terms corresponding to arbitrary 
'subtraction constants'. 

If composite operators need to be defined, they may require new coun- 
terterms in order to recover gauge invariance. See section 6 for the example 
of global currents. 

Finally, amplitudes which involve external fermion lines can continue to 
be treated using the x w propagator in the background of the gauge field. 

3 Gauge Field Interpolation 

We describe our adaptation of a recent suggestion by 't Hooft ||] for inter- 
polating a gauge field living on the links of the 6-lattice. We consider the 
two cases of interpolation to any /-lattice subdivision of the 6-lattice, or to 
the continuum. The latter will be of use since many of the limiting proper- 
ties of the /-lattice interpolation as / — > will be those of the continuum 
interpolation. After having defined the interpolation algorithm we discuss 
the important properties that it satisfies. An alternative interpolation which 
may be computationally simpler is discussed in section 8. 

Broadly speaking, the interpolation is obtained by pasting together so- 
lutions of the continuum or /-lattice Yang-Mills equations in each four- 
dimensional hypercube of the 6-lattice, where the gauge field on the bound- 
aries of the hypercube are specified in the directions tangential to the bound- 
ary. This should determine a unique gauge field in the hypercube up to gauge 
transformations in its interior, which we fix with a Lorentz gauge condi- 
tion]]] The tangential gauge fields on the three-dimensional boundary cubes 

^The existence of solutions and their uniqueness once the Lorentz gauge condition is 
imposed can be proven to all orders in the gauge field self-coupling. This result should 
apply when the flux through any 6-plaquette is sufficiently small, as will be the case near 
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are themselves determined by three-dimensional Yang-Mills equations, and 
similarly for their own boundaries, the two-dimensional plaquettes. The 
continuum or /-lattice gauge fields tangential to the links of the 6-lattice 
are fixed to be constant along the link, and these set the ultimate bound- 
ary conditions for the multi-stage interpolation procedure. The details now 
follow. 

Let U^(s) = exp(zL4 M (s)) be the gauge field configuration on the 6- 
lattice. We can describe the interpolation to either the continuum or the 
/-lattice in terms of a field a«(x), which represents the continuum gauge 
field, or the /-lattice gauge field via u^(x) = exp(z/a M (x)) (in which case x 
is taken to be restricted to the /-lattice vertices). First, on points lying on 
the links of the 6-lattice we interpolate the tangential gauge field component 
by, 

a ll (x = s + b\fi) = A ll (s), (24) 

where < A < 1. The interpolation into a 1-2-directed 6-lattice plaquette is 
given by minimizing the two-dimensional Yang-Mills action 

S = J d Xl dx 2 ff 2 (x), (25) 

or 

S = J2 Tr [! " Mx)M* + n>i\x + fl + f2)u^(x + /2)] + h.c, (26) 

X 

in the gauge d\a\ + d 2 a 2 = 0, subject to the boundary conditions given 
by agreement with the a\ and a 2 already determined on the bordering 1- 
and 2-directed links respectively. The x being integrated (summed) over are 
restricted to the points inside the 6-lattice plaquette. f^ v denotes the field 
strength made from a^. Note that only the a M components tangential to 
the 6-plaquette have been determined. The generalization to interpolating 
in 6-plaquettes pointing in other directions is obvious. 

The way to complete the interpolation procedure is also clear. We in- 
terpolate to the interior of any 6-lattice 3-cube by minimizing the appropri- 
ate three-dimensional Yang-Mills action in three-dimensional Lorentz gauge 
with the boundary condition being agreement with the tangential gauge 
field on the plaquettes bordering the cube, determined at the previous stage. 
This will determine the three components of the gauge field tangential to the 
cube. Finally we fully interpolate the gauge field to every four-dimensional 

the continuum limit. 
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hypercube by minimizing the four-dimensional Yang-Mills action in Lorentz 
gauge with the boundary conditions being agreement with the tangential 
gauge fields in the 3-cubes bordering the hypercube. This completes the 
determination of a^(x). 

This interpolation procedure is a finite computation (in finite volume) 
when going from a 6- to /-lattice. It satisfies several important properties 
which we now list and briefly discuss. 

• Transverse continuity: This describes how smooth our interpolation is. 
As one passes from one hypercube of the 6-lattice into another, the compo- 
nents transverse to the direction taken (ie. tangential to the 3-dimensional 
boundary) vary continuously because these components were fixed on the 
boundary as part of the boundary conditions for interpolation into each 
of the adjacent hypercubes. However the 'longitudinal' component of the 
gauge field can vary discontinuously. This is also true at lower dimensional 
levels: passing between adjacent 3-cubes or adjacent plaquettes or adjacent 
links of the 6-lattice. 

One may wonder what the relevance of transverse continuity is. The 
point is that the interpolation procedure does not sharply cut off the mo- 
menta of the fields above n/b even though of course the 6-lattice fields 
we are interpolating have sharply cut off momenta. Rather, the interpola- 
tion procedure should be seen as a 'soft' cutoff of the a M momenta above 
roughly 1/6. The less smooth the interpolation, the less sharp the cutoff of 
a M momentum. Now suppose T[a] is a well-defined functional of continuum 
gauge fields a^. We will want r[a[£/]] to be a finite functional of U. This will 
fail unless a[U] momentum components fall off fast enough in the ultravio- 
let. (In an /-lattice we will not get an infinity but will get unsubtractable 
dependences on 1// which will preclude our taking the / — > limit.) This 
important consideration is taken up in detail in the next section where it 
is shown that transverse continuity is sufficient to ensure a well-behaved 
fermion determinant as a functional of U^. 

• Gauge Invariance: A gauge invariant functional of the interpolated 
gauge field, T[a], automatically yields a gauge invariant functional of the 
6-lattice gauge field, r[a[f/]. Clearly, it is sufficient to show that 

a[U n ]=a u [U], (27) 

where O(s) is a 6-lattice gauge transformation of U^, and uj[U,Q](x) is a 
continuum (/-lattice) gauge transformation which interpolates f2 and can 
depend on U. We show that such an uj exists by describing its construction 
in the links, then plaquettes, cubes and hypercubes spanned by the 6-lattice. 
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For points on a link of the 6-lattice in the [x direction we take 



u{ x = s + A/26) = (^(s)C/ At (s)0" 1 (s + fib)r x n{s)U ll (s) x . (28) 

This has the property that a^(s + A/26)) = A n (s), where = e lA ^ , which 
is indeed the behavior of a[U ] on points on the b- link. Inside a 1-2 directed 
6-lattice plaquette we extend uj defined on the bounding links in a continuous 
fashion, so that a? is now defined in the plaquette. Since one is always free 
to take a gauge field to Lorentz gauge via a gauge transformation, uj can be 
chosen so that a u is in the two-dimensional Lorentz gauge d\af + = 0. 
By the gauge invariance of the two-dimensional Yang-Mills action, a" yields 
the minimum action of gauge fields agreeing with the tangential components 
a w = AF 1 already defined on the bounding links. These are precisely the 
defining properties of the interpolation a[f/ n ] on the plaquette, so a w = 
a[U ] there. Again, the remainder of the construction should be clear; uj can 
be extended into cubes and hypercubes till one has a gauge transformation 
everywhere and a? satisfies the defining properties of a[U n ]. 

• Locality: When the momentum space Feynman vertex for the coupling 
of the interpolation a p to the fermions is expressed in terms of An, the result 
is a series of interaction vertices which are polynomial in the A p momenta 
and fields, with higher dimension interactions balanced by powers of b. For 
low momenta q <C 1/6, 

a^q) ~ A^q). (29) 

This allows us to view the whole interpolation procedure as an intermediate 
step in coupling 6-lattice gauge fields to fermions that live in the continuum 
(/-lattice) in such a way that the low-energy physics far below 1/6 is unaf- 
fected. The complexity of the procedure arises from the need to maintain 
gauge invariance. 

The reason locality is satisfied by the interpolation is quite easy to see: a p 
at the point x = s + yb is determined purely by A fl (s),A fl (s + ub), that is by 
A^s), d u A^(s), and agrees with A p on the original links of the 6-lattice. For 



low momentum this corresponds to eq. ( 2£ ) . Overall momentum conservation 



is violated by multiples of 27r/6 as is familiar in lattice theories in general, 
the difference in the present context being that the fermion momenta live 
in the /-lattice Brillouin zone not the 6-lattice Brillouin zone. This is not a 
problem for locality but is related to taking the / — > limit, which is dealt 
with in the next section. 

• Rotational and translational symmetry: Our interpolation scheme is 
manifestly invariant under rotations and translations of the 6-lattice so any 
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invariant functional of translates into an invariant functional of 

• Topological Winding Number: As noted by 't Hooft ||, this type of 
interpolation allows us to ascribe a winding number to any 6-lattice config- 
uration whose field strength dies off at infinity, namely the winding number 
of its continuum interpolation. This will be useful in section 7, dealing with 
fermion number violation. Clearly, 'instantons' smaller than b cannot occur 
as the interpolation of a 6-lattice gauge field. The basic idea of assigning 



winding number to lattice gauge fields was first explored in ref. [18|. 

It is good to see the simplest concrete example of many of the above 
considerations. This is provided by the interpolation of U(l) gauge fields, 
for which we can explicitely construct the solution, 

ai (x = s + yb)= Ms + b^2z p p) l[(z„(2y v -l) + l-y u ), (30) 

22,23,24=0,1 p^l v^X 

the other components obtained by cyclic permutation. It is easily checked 
that this interpolation satisfies all the right boundary conditions, is in Lorentz 
gauge and satisfies Maxwell's equations (minimizes the Maxwell action for 
the appropriate dimension) in the interiors of hypercubes, 3-cubes and pla- 
quettes, and is transversely but not longitudinally continuous. Also, Mattice 
gauge fields which are gauge equivalent interpolate to gauge equivalent fields 
in the continuum or /-lattice. 

It is also instructive to examine the interpolation in momentum space, 

where we decompose continuum (/-lattice) momenta according to 

27T 

Qn = Qft + ~b m v- wnere l<?/J < n /b, and is an integer. (32) 

We see that a M has momentum components above ir/b (m„ 7^ 0), though 
does not. But because of the transverse continuity of the interpola- 
tion, above ir/b the momentum components fall off (with m„), faster in the 
transverse directions than in the longitudinal direction. This type of falloff 
will be sufficient for our purposes in the next section. At low momentum, 
q u b <C 1, the interpolation looks very simple, a^(q) ~ A^q). Finally, note 
that the interpolation can be expanded as a power series in bq^ so that there 
is no non-locality introduced by substituting the interpolation into a local 
interaction vertex of a^{q) to fermion fields. 
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4 The / Limit 



We have seen in section 2 that r r [u] at one loop is finite in the limit / — > 0, 
and gauge invariant up to 0(f 2 ). Also, it was shown that G[u] is gauge 
covariant up to 0(f 2 ). However, when T r and G are considered as function- 
als of the real dynamical variables, U, since the interpolation u[U] depends 
on /, we must take care that no extra powers of 1// are generated. The 
objective of this section is then to prove, 

limf-,o T r [u[U]] < oo, lirrif^o G[u[U\] < oo, (33) 

and 

r r [ u [u n }}-T r [ u [u}} = o(f 2 ), 

G^ 1 HU n ]}-G[u[U}}=0(f 2 ), (34) 

where Q(s) is a gauge transformation on the 6-lattice, while uj[U, 0] is its 
interpolation to the /-lattice. ^Recall that our interpolation has the property 
that = u w [U].) G^ is the obvious gauge transformed fermion 

propagator. 



4.1 Fermionic Effective Action 

In perturbation theory, the renormalized one loop effective action in terms 
of the interpolated fields is given by, 

r >] =E / I[d%r^ a2 j qi ,.., qn ;f) aM-^aM <5WQ>»(35) 

n JBZ U) i i 

which can be rewritten as, 

T r = J2 / U d ^i E rW 2 ...(gi,..,gn;/)a Ql ( ( Zi)...a an ( (?n )5( 4 )(^ % )(36) 

where qi = q~i + 2irmi/b and m w = -b/f, .., b/f. 

Although the amplitudes are finite in the limit / — > 0, it is clear 
from eq. (|36|) that the sums over integers m w could introduce powers of 
1//. Obviously, for a smooth enough interpolation, the high momentum 
components of the interpolated fields are expected to be cut off by ~ 1/6, 
in such a way that the sums over m w are convergent in the limit / — > 0. 
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We now prove that this is indeed the case. For the proof we only need 
some smoothness properties of the interpolation and a bound on one loop 
(fermion) lattice integrals for large external momenta q ~ 1//. 

In the previous section we constructed a tranversely continuous interpo- 
lation. The gauge field however jumps in the longitudinal direction at the 
boundaries of the 6-cell, i.e. a fl (n) is not continuous in the fi direction, but 
it is in any tranverse direction v ^ a. Also, within each 6-cell it is fully 
continuous with continuous derivatives to arbitrary order. By continuity on 
the lattice we mean the following property, 

a M (n)-a>-l) = 0(/). (37) 

Note that this definition already implies the existence of a derivative that 
maybe discontinuous, but is bounded for small /. For an interpolated field 
with these properties, it is easy to derive a bound on the high momentum 
components, 

M*)l < C]tt II lir < ^itt II i ( 38 ) 

with 

% = j • (39) 

where C, C are constants independent of /. Note that this bound is satisfied 
by the U(l) interpolation of the previous section. 

We sketch the proof for a function of one variable, the extension be- 
ing straightforward. A function, F(n), defined on the points nf, which is 
not continuous at the boundaries nf = sb (but continuous with continuous 
derivatives elsewhere) can always be written as the sum of a step-like func- 
tion S (i.e. constant in each interval (s — 1)6 < nf < sb) plus a continuous 
function F, 

F(n)=F(n) + S(sj), (40) 

with nf = sb + Snf, and < Sn < bj f . The momentum components of 
S behave as 0(l/|g|) < C/\q\, for q ~ 1/f. On the other hand, since F is 
continuous in the sense of (|37|), its first lattice derivative is bounded in the 
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limit / — > 0. Also the momentum components of this derivative are bounded 
and we can write, 

1^)1 = ,0«5M| < c > (4i) 

This implies that |-F(g)| < C/|<z|, with C independent of /. It is easy to 
repeat the analysis for a continuous function, by writing its derivative as 
the sum of a step-like function and a continuous function. We find in this 
case \F(q)\ < C/\q\ 2 . 

The second ingredient in the proof of finiteness is a simplified version of 
the power counting theorem of Reisz . We use the derivation of reference 
[ppf for one loop diagrams. Here, we simply state the result and refer to 
appendix A for the details. The renormalized lattice integrals contributing 
to the r(") of @ satisfy the bound 

\T^(q;f)\ < C\q\ 4 - n Logs(\q\), (42) 

for all f < fo for some fixed /o, with C a constant independent of /. n is the 
number of external boson lines and q denotes generically the external mo- 
menta. Note that this is also the result expected from naive power counting 
in the continuum. 



Now, we are ready to prove the finiteness of (36). First we note that 
for n > 4, eqs. (|38| ) and (|4^) imply that the sums over are absolutely 
convergent in the limit / — ► 0. Since the genuine logarithmic divergence in 
the n = 4 amplitude has been absorbed into the definition of gj, (see eq. (|l8|) ), 
this means that all the contributions with n > 4 are finite in this limit. 

To prove finiteness in the n = 2 and n = 3 cases, we have to rely on the 
transverse properties of these amplitudes for the hybrid action. The P-even 
contributions are gauge invariant. This implies, 

£&rg>(9) = ° EC 2 ri t lfe,9 2 ) = cr( 2 c )fe il ). (43) 

It is shown in Appendix A that these relations are enough to prove the 
stricter bounds, 

|r$(?)l < C\q a \ 2 Log(\q\) \r^(q)\ < C\q a \Log(\q\) a + M . (44) 



Together with (pq ) and (|42D , it is now straightforward to prove absolute 
convergence of the sums for the parity-even n = 2 and n = 3. 
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On the other hand, these relations are not satisfied by the x w action 
(0) and in fact we cannot prove convergence of the sums in (36). This 



is simple to understand in real space. We have seen that the chiral Wilson 
action requires a counterterm of the type, 

S X (f) Tr ]T(A M u M ) 2 . (45) 

This is not well defined in the / — > limit, because the longitudinal deriva- 
tive of the interpolated field is a sum of approximate delta functions, whose 
square is ill-defined. For finite /, it would diverge as a power of b/f, intro- 
ducing more divergences than those in Z\y. In the absence of a smoother 
interpolation, we conclude that the chiral Wilson action (]?]) cannot be used. 

Finally, we need to consider the parity-odd parts in n = 2, 3. The P-odd 
amplitude for n = 2 is exactly zero. On the other hand, the three-point 
function is not zero. Obviously, since anomalies cancel, these contributions 
are gauge invariant in the limit /—►(), but for finite /, we cannot use 
the fact that the diagram is gauge invariant to derive a relation like fl44}). 
Nevertheless, this bound still holds because of the presence of an epsilon 
tensor. It is proven in appendix A by direct inspection. Again this is enough 
to prove absolute convergence of the sums. 

This concludes the proof that the one-loop renormalized fermionic effec- 
tive action derived for the hybrid scheme (^) has a good limit / —* 0. 

We still need to prove that the gauge non-invariant contributions that 
remain after renormalization are 0(/ 2 ). In section 2, we have seen that this 
is true for r r [u]. Here we prove that, when the / dependence of the interpo- 
lation is taken into account in rV [«[£/]], no extra powers of 1// change this 
result (U|). For this we just need to extend the bounds ( [42|) to amplitudes 
with one insertion of 5 U S^ W . We denote these amplitudes by d^T^. Under 
a gauge transformation on the 6-lattice, r r for the hybrid scheme changes 
by, 

det[iD\* J 

r b/f 

E / EK* E ^FW( gi )afe)afe)...a(g„)5 (4) (E*)' ( 46 ) 
n jBZ ( b ) i m, H =-b/f i 

where, S* w is given by the action (0). In momentum space, we get 

S u gf = i9 a (q)5 {4 \q + p-p')y(p'){M(p')L - M(p)R}T a ^(p), (47) 
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for an infinitesimal gauge transformation uj = e t9a ( n ) T °.^ obtained by inter- 
polating an infinitesimal ^-lattice gauge transformation f2(s) (see section 3), 
and M(p) = (2r//) J2 V s i n2 (py//2). Inserting this term into a free fermion 
propagator gives, 

M(p>) 2 1 M{pf 1 

H pj{ P ')p f (v) Pj(p)P f (pT ( j 

where 1/-P/ is the free fermion propagator. Also, we have the bound, 

|^M|<Cy (49) 

for all f < fo and C some constant. From the analysis of Appendix A, we 
conclude that a diagram with one insertion of 5 w Sf and n gauge legs, is the 
sum of local terms up to dimension 5 plus a non-local part which is bounded 
by, 

\5^\<f 2 \qf- n Logs(\q\). (50) 

Since there is no local P-odd operator of dimension 5 or less, except the 
anomaly, which vanishes, the non-local part bounded by ([5(]) is all there is 
and as expected is 0(f 2 ). 

It is clear from the previous discussion that amplitudes with n > 6 have 
two or more powers of f 2 and since for large momentum, q ~ 1//, they 
must satisfy boun d (|50[ ), we can use it together with (|38|), to prove that the 
sums over in (jig ) are absolutely convergent and do not bring any extra 
powers of 1//. 

For n < 5, we make use of the fact that there is an epsilon tensor. In 
the case of n = 2, it can be shown by direct inspection that the integral is 
bounded by, 

KT ( - 2 Xq)\<Cf 2 \q a \\q^\\q\ 2 Log(\q\) a,P^fi,u a + (3 (51) 

Using (^) and the full continuity of the interpolated gauge transformation 
0(n), which leads to a bound of the form 

< C U^^ (52) 

a "a 

it is easy to see that the sums over yield a result suppressed by f 2 up 
to logs. A similar analysis shows this is true for n = 3 — 5 diagrams. This 
ends the proof of (pi) . 
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4.2 Fermion Propagator 



We finally need to consider the propagator G [it [[/]]. We will prove that it 
has a well-defined / — > limit and that the breakage of gauge invariance is 
down by at least 0(f 2 ). In perturbation theory, it is given by, 

/n n 
II dpi J2 G {n) (Pi, ..,Pn)a{px)...a(j> n )5^{q' - q - 5»(53) 
« 3Z ( b ) i=i m Mi i=l 

where as usual Pi = pi + 2-Kvrii/b. Since, the amplitudes are products 
of free propagators, and \Pf(k)\ < C/\k\, this implies that for large m^, we 
must have |G?( n )| < C/\pi\, for any large pj. This bound together with (|38|) 



implies the absolute convergence of the sums over in (53). The / — > 
limit of G[U] is well-defined. 

Now we have to prove that the gauge non-invariant contributions to these 
diagrams are 0(f 2 ). The gauge variation of G[U] is given perturbatively by 
a sum of diagrams with an arbitray number of gauge boson legs and one 
insertion of 6uSj w . Using ( |47| ) and (fl"9|), we see that the line with the 
insertion can be bounded for large by f 2 \k\ where k is the momentum of 
the fermion line, as compared to the behaviour ~ 1/|&| without the insertion. 
From (|53|), it is clear that the sums over are at most logarithmically 
divergent in this case, so that the breakeage of gauge invariance in G[U] is 
at most 0(f 2 Log(f)), as expected. 



5 Gauge Ward Identities. 

We are now in a position to demonstrate the gauge invariance of our con- 
struction when gauge fields, U, are integrated in the functional integral. 
We have shown in the previous section that the three components of the 
full functional integral (see eq. @), T r [u[U]],G[u[U]],S g [U], satisfy gauge 
invariance up to 0(f 2 ) (in fact we showed that mathematically the gauge 
invariant limit / — > exists, but in any simulation / will be small and non- 
zero). In the absence of gauge fixing, the full generating functional satisfies, 

"* s 

I WdU^s) e -s.M+£.^, e ?rHU]] e -*J[«Mto + 0( A (54) 

"'a 
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where is any 6-lattice gauge transformation while uj[Q, U] is its /-lattice 
interpolation described in section 3. The gauge transformed propagator is 
given by, 

G^HC/^Kn!,?^) = u- 1 (n 1 )G[u[U n }]{m,n 2 )u(n 2 ). (55) 

The central observation needed to understand eq.(|5"4|) is that since the U 
fields live on the 6-lattice, the largest momentum scale they can introduce 
when integrated is 0(1/6). That is, they cannot reintroduce any powers 
of 1/ f which could cancel the f 2 suppression of gauge symmetry violation. 
Note that this is a non-perturbative expectation |JJ| 

For perturbative investigations it is useful to add the standard, rotation- 
ally invariant gauge-fixing and ghost terms, neglecting the small breakage 
of gauge invariance in r r [it[t/]] and G[u[C/]] in the Fadeev-Popov procedure. 
This results in lattice BRST identities among the Green functions which are 
broken only by 0(/ 2 /6 2 ), since the gauge boson loop momenta are cut off 
above 0(1/6). The statement of the gauge invariance and unitarity of the 
quantum theory is contained in the BRST identities. From these it should 
also be possible to prove perturbative renormalizability and the continuum 
(renormalized) BRST identities of the theory. Of course, it is believed that 
there is a unique continuum theory with the particle content of our lattice 
construction which satisfies the BRST identities pertubatively, and non- 
perturbatively coincides with the non-gauge-fixed theory. 

Note that to compute physical amplitudes, it is essential to be able to 
form gauge invariant combinations of the external fermion fields and the 
gauge fields. This, at first sight, seems problematic since the fields live on 
different lattices, but clearly the solution is to make use of the gauge field 
interpolation u[U] in forming path-ordered exponentials of the gauge field 
connecting external fermions. 

Let us turn to the question of how small to take the ratio //6 in a 
lattice simulation in a (physical) volume, L 4 . Clearly we should take f/b 
small enough that the measurements in our volume are insensitive to the 
breakage of gauge invariance. For example, a simple perturbative estimate 
suggests that (in the gauge-fixed theory) the 0(/ 2 /6 2 ) breaking of the Ward 
identities can induce a gauge boson mass of order //6 2 . Clearly, if//6 < b/L 

'Small deviations from exact gauge invariance, such as a small bare gauge boson mass 
term, lead to difficulties in continuum perturbation theory, as is well known. Indeed one 
is really outside the perturbative regime. However non-perturbatively there is no sign of 
difficulty. For a brief review of the lattice literature on the example of bare gauge boson 
masses see ref. [pit] . 
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we would be insensitive to this violation of gauge invariance in our volume, 
since the gauge boson correlation length would be larger than L. With such 
a relationship of scales, in the continuum or thermodynamic limits, b/L — > 0, 
gauge invariance becomes exact (/ jb — * 0). 



6 Anomalies in Gauge and Global Currents 

In this section we will first check that our hybrid scheme reproduces the 
consistent gauge anomaly at the one-loop level. As explained in section 
2, this is to be expected on general grounds. Obviously, we only consider 
theories where the gauge anomaly vanishes (i.e. d a h c among the fermion 
representations is zero). However, in most theories of interest, there are also 
other currents associated with global symmetries that can be anomalous 
(e.g. B+L number in the Standard Model). We will show that these are 
also properly reproduced. In studying current conservation, we will need to 
consider Green functions with one insertion of the global current. In general, 
new counterterms are needed to ensure gauge invariance of amplitudes con- 
taining composite operators such as this current, even in the hybrid scheme. 
Once these counterterms have been added, we show that our construction 
gives the correct gauge invariant anomaly at the one-loop level. It is then 



straightforward to see that the Adler-Bardeen theorem [22] must hold to all 
orders in perturbation theory. 

The axial anomaly triangle in a vector-like theory was computed in lat- 
tice perturbation theory with standard Wilson fermions by several authors 
Hi 23]. The case we have to consider here is slightly different, since we are in 



a chiral gauge theory q The gauge anomaly triangle for our hybrid scheme 
is given by the imaginary part of the diagram with one insertion of ( [47] ) and 
two gauge boson legs, in the limit / — > (recall that in the hybrid scheme 
only Im V breaks the gauge symmetry). The result we get is, 



5uF$ bC {q,q) = -^0 a (q + q')d abc e^ aL ,f3q' a q{3 x 

4 2 rMs 2 -cpM 2 2 

d k c^c a c uS/3 ^ 2 ~ M 2) 4 + 0(f ) (56) 

where w = exp{i9 a T a ). Also c a = cos(k a ), s a = sin(k a ), M = 2rJ2 a sin(k a /2y 
and s 2 = Xla s ™(^a) 2 ' r is the Wilson parameter. This integral can be 

**For a discussion of gauge and global anomalies in the context of the 'overlap' formu- 
lation of chiral gauge theories see ref. VM 
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solved exactly Q and gives a result independent of r as expected, 

6^ 2 \q,q>)= -^e a (q + q>)d ab % auP q' aqi3 + 0(f 2 ). (57) 



The Wess-Zumino consistency conditions [12] then give the second term of 
the gauge anomaly, 

<Lr = -j2-^e a M3vTr[6 d a (a^dpa u - -a^apa,,)] + 0(/ 2 ). (58) 

We have proven earlier that, in the absence of gauge anomalies, the hy- 
brid scheme gives a gauge symmetric continuum limit. However, the chiral 
gauge theories we want to simulate typically contain global t/(l)'s which 
are anomalous. We will consider a general case, in which a global current 
is denoted by jn, with generator t. We will assume that t commutes with 
the T a . A convenient way to find the anomaly is introducing an external 
source b^(x) that couples to the global current. It also lives on the links of 
the /-lattice: = exp(ifb^t). The simplest choice is to couple it like the 
gauge field (any other choice would lead to the same result). In the hybrid 
scheme we just substitute, 

Uf, ->• u^u^ (59) 

both in the chiral Wilson action (0) and in the standard Wilson action, 
defining det[iD] xw and det[i p\ sw in (|23|), 

Note that for the chiral Wilson action this prescription corresponds to 
only the physical left-handed fermion field transforming under the global 
symmetry. 

Due to the presence of the external field u, the fermionic effective action 
is no longer gauge invariant for d a bc = 0. Instead we get, 

M>,«] = -r^e a ^ v Tr{6d a {2 a„d p b v - ^a p a u )] + O(ftJ) + 0(/ 2 ) (60) 

Since we are only interested in one insertion of the current, we just need the 
linear terms in b^. It is easy to see that by adding the local counterterm, 

i 3 
S c .t. [a, b] = j^-^e aflf 3 u Tr[2 b a a^dj3a u - i-baa^a^] (61) 



^We thank Ph.Boucaud for pointing this out to us. 
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one recovers gauge invar iance, 



Mr + s c .O = o(f). (62) 

Finally, after this subtraction, the global anomaly is given by, 

8\(T + S c .t.) \b^=o = Im{< SxSf > xw )\ b/i=0 + 6\S c . t . \b„=o (63) 

where A = e 1 ^^ 1 . The first term corresponds to a diagram which is identical 
to the one for the gauge anomaly, with </> t replacing 6 a T a in (58), 

S X T | 6=0 = -^2^2^ ea ^ uTr ^ dotia^dpCLv - -a^apa u )} + 0(f 2 ) (64) 
When we add the counterterm we find, 

4> = 5 X (T + S ct ) \ h=0 = -J^2^ Tr[t f^f^] + O(f) (65) 

with f^y being the gauge field strength. This is the expected gauge invariant 
form of the anomaly. 

This is the result at one loop. However, since gauge loops can at most 
generate extra powers of 1/6, all the 0(f 2 ) effects in ( |65| ) go to zero in 
the continuum limit, f /b — > 0, to all orders in perturbation theory. This 
means that in the continuum limit, we recover the well-known result of the 
Adler-Bardeen theorem since, 

<d lx j ll F> = < --l-^Tr[tf^f fll/ }F> + <5 t F>+0(f 2 /b 2 ) (66) 

where F[U, is an arbitrary functional. 



7 Fermion Number Violating Amplitudes. 

Consider the example of a theory with one generation of standard model 
fermions interacting with gauge fields but no scalars (since we have not 
considered them in this paper). In this section we will understand the 
phenomenon of B + L violation j2f| in an 'instanton gas' picture, in the 
continuum and then on the lattice. This will be sufficient to demonstrate 
that B + L-violating processes will occur in our lattice construction. We 
then give a prescription for defining lattice-regulated net-B + L violating 
amplitudes. A similar discussion has been previously given in reference |24fl , 
from the point of view of the 'overlap' formulation of chiral gauge theories. 
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The main problem is that we have restricted the lattice construction to 
net fermion number conserving amplitudes (for every external *$>l we have a 
corresponding external ^>l)- Still, as we saw in the last section, even in this 
sector of the theory there is local B + L violation, as given by the anomaly. 
This will allow us to infer B + L violation in the full theory. The strategy for 
defining more general correlators, say < 0{x) > representing the amplitude 
for a B + L- violating process in the vicinity of x, is to demand cluster 
decompostion of the full theory and consider instead the fermion number 
conserving correlator [|l0| , 

<0{x)0\y)>\ x ^ O0 ~ <0(x)><0\y)>, (67) 

which is defined in our construction. That is, if this two-point correlator 
does not vanish as x and y are taken arbitrarily far apart, then B + L is 
being locally violated (with opposite signs) at x and y. Since events at y are 
too distant to affect those at x in a sensible theory, we must have that 

< 0{x) 0. (68) 

Inferring fermion number violation by considering fermion number conserv- 
ing correlators has been exploited in refs. p6|]. 

We will first see how this works formally in the continuum, using the 
formalism of Fujikawa for chiral theories p7[ |, and then consider our lattice 
construction afterwards. For simplicity, we will take 

Oix) = q L {x)q L (x)q L (x)£ L (x), (69) 

where we take the internal indices to be contracted gauge invariantly. This 
operator carries two units of B + L. We then have 

oHy) = q L (y)qL(y)q L (y)My)- (70) 

For a fixed background gauge field a^, we can decompose the chiral quantum 
fields in terms of the complete, orthonormal set of eigenfunctions of the 
vector-like hermitian Dirac operator, i I/), 

i ]fiip n = A n VVi 

n>0 

^ L = E««^ f - ( 71 ) 

n>0 
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Here, we are using a condensed notation where really there are separate 
Dirac equations for the = Ql,^l fields. The two-point correlator of O is 
then given by functional integration over the ferniion coefficients Q>ni 

< 0(x)0\y) > a = (I] KfiU AlK Imr (£ g " (x y ( ^ ) 3 (E ^ (X) 5 t(y) : 

n>0 n>0 n>0 " n>0 n 

(72) 

where < ... > Qm denotes the fermionic averaging done in the background 
of a„, which is presently unintegrated. The origin of the various terms is 
as follows. The products of eigenvalues for the various fields comes from 
integrating the fermion action, which is diagonalized in this basis. This is 
nothing but the square root of the vector-like determinant, which we gauge- 
invariantly regulated in section 2. The phase factor emerges as a Jacobian 
for a formally unitary transformation, 

f V^ L V^ L = / II ^ n da n e iImr . (73) 

J J n>0 

The bilinears in the eigenfunctions are the fermion propagators in the eigen- 
basis. 

Now we specialize to the case of a zero winding number background 
gauge field which is made out of a distantly separated instanton and anti- 
instanton pair, + o^, the instanton in the vicinity of x and the anti- 
instanton in the vicinity of y. To understand the dominant contribution 
to our two-point correlator we note that for a single instanton, the Atiyah- 



Singer index theorem [28| tells us that each of the fields comprising O has 
a left-handed zero mode of its associated vector-like Dirac eigen-equation. 
For an anti-instanton background, there are right-handed zero modes of 
the Dirac eigen-equation instead. In the background of the instanton anti- 
instanton pair, the left and right-handed zero modes pair up to form a Dirac 
spinor eigenfunction with non-zero eigenvalue, Ai. As the instanton and 
anti-instanton are moved far apart, Ai — > and the two chiralities separate 
spatially. As this happens the fermion propagators are dominated by the Ai 
modes, 

h — ^ ^ — ' ( } 

n>0 " 1 



while the product of eigenvalues is suppressed by Ai- Cancelling the powers 
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of Ai gives 

< 0{x)0\y) > a ^ ~ (I] A^) 3 (I] X e n )e^ r (qi:(x)qf\y)f(ef(x)lf\y)), 

n>l n>l 

(75) 

where the leading behaviors of ty\ and ipf^ are given by the zero modes of 
the instanton and anti-instanton respectively. 

The mechanism by which < OO^ > survives as \x — y\ — > oo, when we 
integrate over gauge fields in a dilute instanton gas approximation, is now 
clear. There will be a purely perturbative contribution to the correlator 
which will fall to zero as \x — y\ — > oo and a contribution from an instanton 
anti-instanton pair of the above form which does not fall off, but rather 
factorizes into contributions from the vicinities of x and y. We are forced to 
conclude that B + L violation occurs in the full theory. Of course, we expect 
the qualitative point we are making to survive the (rather drastic) dilute 
gas approximation. The moral is that the physics of non-trivial winding 
number is implicit in the zero winding number sector and this is linked to 
B + L violation through the (integral of the) anomaly equation for the B + L 
current. 

Let us now see how all this relates to our lattice construction. As de- 
scribed in section 3 we mean by a single 'instanton', a 6-lattice configuration 
whose interpolation to the continuum yields a unit winding number contin- 
uum gauge field. This is easily arranged by obtaining the original 6-lattice 
gauge field by discretizing a large (compared with b) continuum instanton 
in the first place. An instanton and anti-instanton are then easily composed 
to give a zero winding number gauge field. For / <C b <C instanton size, the 
hybrid regulator gives essentially the same result as our continuum anal- 
ysis except that the product of eigenvalues is cut off above roughly 1//. 
The vector-like determinant is determined by the gauge invariant standard 
Wilson action. We can work with the related set of lattice eigenfunctions, 
which for small / we expect to approximate the Fujikawa eigenfunctions of 
our continuum discussion. The gauge field interpolation procedure makes 
the remaining analysis as simple as the formal continuum case. The fermion 
propagators are nearly the same as the continuum propagators in the back- 
ground of the interpolated gauge field. The poles corresponding to the 
doubler modes have completely decoupled. Therefore the O two-point cor- 
relator is dominated by the approximate zero modes just as in the earlier 
analysis. For \x — y\ —> oo the correlator must be interpreted as B + L vio- 
lation occurring at x and at y but without any physical connection between 
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the two points. 

The above discussion leads us to a straightforward means of defining 
arbitrary lattice-regulated fermion number violating amplitudes in terms of 
the fermion conserving amplitudes we have already defined. For example, 
let O be a product of fields at various points, which carries two units of 
B + L. Then for x far removed from all these points, cluster decomposition 
yields 

< 6o\x) > ~ < 6 >< 0\x) > . (76) 
We can satisfy this by defining 

< O > = — r , (77) 

(< 0{x)&(y) >)s 

for y, x also distantly separated. Note that all the quantities on the right- 
hand side are defined in the zero- winding number sector^. 



8 Computational Efficiency. 

Obviously, the most computationally expensive operation in a simulation of 
chiral gauge theory using the methods of this paper, compared with lattice 
QCD say, is the calculation of the fermion determinant on the /-lattice, 
instead of the coarser 6-lattice. The other major complication is the gauge 
field interpolation technique. Let us see what can be done to reduce these 
costs. 

In the hybrid scheme described in section 2, there are in fact two de- 
terminants, the sw and x w ones, which must be calculated on the /-lattice 
for each gauge field configuration. The consistent use of the /-lattice for 
fermions was done just to reduce confusion. It is however possible, and much 
more efficient to calculate the sw determinant directly on the 6-lattice! The 
point is that the sw determinant is exactly gauge invariant and we know 
how to non-perturbatively restore the global chiral symmetry broken by the 
standard Wilson term on a single lattice for gauge fields and fermions. Just 
as in lattice QCD, we can use the action 

b 3 

s,n 

** Actually, there is also an undetermined phase in < O >, which we have taken as zero. 
Its value in general defines the 6- vacuum we occupy. 
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+ (4r + M6)6 3 ^$(s)$(s) 

s 

- ! ^E[ |, ( s )^( s ) <& ( s + / i ) + ^( s + / i ) c/ A: 1<l> ( s )]' ( 78 ) 

S,/i 

where M has been tuned to the point of chiral symmetry for the vector- 
like fermions by tuning the associated Goldstone bosons to be massless. 
We see that in this formulation, there is no need to do coupling constant 
renormalization for distances from / to b. That is, 

9b = go- (79) 

The x w determinant which is used to regulate Imr must still be regulated 
on the /-lattice because it is not exactly gauge invariant, unlike the sw 
determinant. Thus the workload has been cut in half, but not more. 

Similarly, one can use the 6-lattice to regulate external fermions. If a \ w 
action is used on the 6-lattice, the breakage of gauge invariance would not be 
suppressed. But instead, we can use the standard Wilson action above, of 
course restricting ourselves to calculating with the external fermions carrying 
the gauge charge. 

For conceptual purposes, we have already considered matching the two- 
cutoff theory with a one-cutoff theory. Here we note that such a matching 
would result in our ability to compute in a chiral theory with the same effi- 
ciency as in (unquenched) lattice QCD! The Rome group proposal is similar 
to this one cutoff theory; they employ the x w regulator scheme with gauge 
fixing terms and demand that BRST relations among the Green functions 
are obeyed in any simulation, hoping to thereby tune the counterterm coeffi- 
cients to all orders. However, non-perturbatively the viability of this proce- 
dure is unclear. On the other hand, one can imagine computing a handful of 
Green functions using our two-cutoff method, possibly in a volume not too 
much larger than 6 4 . It may then be possible to tune the counterterms in a 
single 6-cutoff theory in order to reproduce the same results to good preci- 
sion. Once the counterterms are determined in this manner they are fixed 
for all other possible Green functions. Therefore one can simply proceed 
to use the single cutoff theory from then on. This corresponds precisely 
to a non-perturbative matching having been performed at b between the 
two-cutoff and one-cutoff theories. 

Next, we observe that the size of the breakage of gauge invariance is 
partly determined by the form of the Wilson term in the x w regulator. The 
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fact that it turned out to be 0(f 2 /b 2 ) is directly linked to the fact that the f- 
lattice Wilson term possessed two lattice derivatives. But nothing prevents 
us from using a four-derivative Wilson term, for example. This will push the 
unsubtractable non-local breakage of gauge invariance to higher order in /. 
If all the lower dimensional local terms are still locally integrable we could 
arrange to have the Ward identities obeyed to higher precision than f 2 /b 2 , so 
a much coarser /-lattice could be used. There is a problem however in that 
the interpolation will almost certainly need to be smoother than it presently 
is. For example, a four-derivative counterterm would be very singular with 
1/f for an interpolation which is only transversely continuous. This kind of 
improvement will have to await the development of a smoother interpolation 
and the check of whether all the local terms that might appear in L(n) are 
locally integrable. 

Finally, the gauge field interpolation as we have defined it is a rather 



cumbersome procedure. A more economical one was given in ref.[29|. It 
satisfies all the criteria we give in section 3, with the possible exception 
of lattice rotational invariance, which we have not been able to prove or 
disprove. However, this symmetry could be restored by averaging over lattice 
rotations, 

r[u] _ Egrp^ (8Q) 

where U R is a lattice rotation of U. It is important to stress that cubic 
invariance is an essential ingredient to recovering a Lorentz invariant con- 
tinuum limit. 



9 Conclusion. 

In this paper, we have presented a non-perturbative lattice formulation of 
chiral gauge theories. It provides a finite algorithm for computer simula- 
tions. The gauge fields live on a lattice with spacing b, taken as the defining 
'cutoff of the theory. For each gauge field configuration the chiral fermion 
determinant is computed by first interpolating the gauge field to a finer 
lattice with spacing / <C b; the determinant is then constructed in the back- 
ground of this interpolated field. Though it is difficult to maintain gauge 
invariance and eliminate the doubler modes that plague lattice fermion con- 
structions, the loss of gauge invariance can be minimized using a separate 
'hybrid' regulator for the magnitude and phase of the fermion determinant. 
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The gauge non-invariance of the determinant is then of order f 2 so long as 
gauge anomalies cancel. The gauge field loops do not upset this suppression 
since their highest momenta are of order l/b. There is one subtlety, that the 
gauge field interpolation is smooth enough to not reintroduce powers of 1// 
that can destroy the approximate gauge invariance. We checked that this 
did not happen with the choice of 'transversely continuous' interpolation 
made here. In the continuum limit / — > 0, followed by b — > 0, exact gauge 
invariance is recovered. 

The price of maintaining gauge invariance is the computational expense 
of gauge field interpolation and calculation of large fermion determinants. 
We argued that in order to be insensitive to the breakage of gauge invariance 
in a simulation in a physical volume L 4 , one would have to take 

f/b ~ b/L. (81) 

Strategies for improved efficiency were discussed including matching the two- 
cutoff theory to a one-cutoff theory. The main computational effort would 
then be in doing the matching. If this is possible, the one-cutoff theory could 
be used for all amplitudes of interest, with the same cost as (unquenched) 
lattice QCD simulations! 

We also explored anomalies in global currents in our lattice construc- 
tion of chiral gauge theory, demonstrating the validity of the Adler-Bardeen 
theorem to all orders. We briefly discussed the related physics of fermion 
number violation, using B + L-violation in a simplified 'standard model' as 
an example. 

At present, the task of using the algorithm presented here to simulate 
four-dimensional, non-abelian chiral gauge theories seems computationally 
formidable. Of course, we hope that there can be improvements made in the 
efficiency of the algorithm, and in the ability of computers to carry them 
out. In the meantime, simulation of two-dimensional chiral gauge theories 
seems quite feasible and provides a useful arena for testing the method. 
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Appendix A 

In this appendix we will prove in more technical detail the general bound 
on one- fermion- loop lattice integrals (p8|), and the stricter bounds (|44|) for 
the gauge invariant parity-even contributions in n = 2, 3 and for the non- 
invariant parity-odd n = 3 amplitudes. 
These integrals have the form, 

r W (q)= [ d 4 kl( g ,k,f)^ [ d 4 k pk'*.;^ n (82) 

where the functions Pf are the lattice fermion propagators, so 

1 



P](k) ^ 



M ( k ) = yE^ 2 fe//2) k^ = -sink^f. (83) 

The fermion momentum is k € BZ(f) = {-j--, j)- 

It has been proven in |2C] for very general conditions on the functions V 
and P, which our action satisfies, that the integrand can be bounded by a 
continuous function that does not depend on /, for all f < fo and /o finite 
and fixed. We define the integer lo, 

V(q,k,f) = f-»V(qf,kf) (84) 

and the degree of divergence of V as, 

limx^ V(q,k\J/X) =K\» + 0{\ V - 1 ) (85) 

and K ^ 0. For the integrals we consider here, the following equality holds, 

2I f -uj = n (86) 

with If the number of internal fermion lines and n being the number of 
external boson lines. The last equality can be easily checked from the lattice 
Feynman rules for fermion loop integrals [O]. 
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With these definitions, Liischer's result is 

\V(q,k,f)\<C\q\^Q"(\k\,\q\), (87) 

where Q is a homogeneous polynomial of degree v and C is a constant 
independent of /. On the other hand, for each propagator we get, 

\Pf(lJ)\ >C'l 2 , (88) 

where C is some constant and I 2 is the square of the fermion momentum. 
Then, 

|r0..(„|</ rt ^f . (89, 

If If is the number of internal fermion lines, we can distinguish two possi- 
bilities. If d uv < 0, 

d uv = 4 + v - 21 f < 4 - n (90) 

the integral is convergent in the limit / — > and by power counting for 
continuum integrals this implies, 

|p(n)| < \ q \^-2I fLogm = \ q \*-n LoM): (91) 

where n is the number of external boson lines. 

In the case d uv > 0, the integral (|89| ) is divergent in the limit / — » 0, 
so it must depend on the cutoff 1//. Note however, that we can separate 
the divergent parts from the finite part by an expansion at zero external 
momenta. From ( |89] ) we see that there are no infrared problems if 

d ir = 21 f - 4 - u < 0, (92) 

but di r = —d U v + v — uj < —d uv . Since the degree of infrared divergence 
increases with each derivative with respect to external momenta, while the 
ultraviolet d uv decreases, we see that for uv divergent diagrams, we can 
Taylor expand up to d uv — 1 times or more. The local terms in the expansion 
are substracted and the remainder is only logarithmically divergent, and also 
possibly non-local. In our hybrid scheme and if the gauge anomalies cancel, 
there is no local counterterm except the one contained in this remainder, 
which is logarithmic and cannot be substracted at zero momentum. Only 
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three cases need to be considered n = 2, 3, 4. In the case of n = 4, d uv < 
so the integral (|89|) is at most logarithmically dependent on the cutoff, 

ir^l <CLo 5 (| g |/). (93) 

We can easily renormalize it by substracting the same integral at zero ex- 
ternal momenta, with an infrared cutoff ~ 1/6, 

<JT (4) (0) = f d 4 kl(0,k,f). (94) 

JBZ(b)->BZ(f) 

(This corresponds to the gauge coupling renormalization of eq. (|lg|).) The 
remainder is then finite. This is easy to prove using the Liischer bounds. 
The renormalized integral is then, 

T^(q) = [ d i kl(q,k,f) + [ d A k{I(q,k,f)-I(Q,kJ)) (95) 

JBZ(b) J BZ(f)—BZ(b) 

The integrand of the second term obviously vanishes for zero external mo- 
menta, so it can be written as a Taylor remainder of a first order expansion 
around zero external momenta. This remainder is of the form, 

R i = E 9a t dt(l - t)^-I(q + tw, k, f)l =0 , (96) 

and now Ri has d uv = —1, so the second term is finite in the limit / — > 0, 
by (p9|). The first term obviously is too. From ( |95"| ) we then find, 

|r( 4 )| <CLog(\q\b), (97) 

which is finite in the limit / — > and has a logarithmic bound for large 
external momenta as expected. 

The cases of n = 2 and n = 3, since they are only logarithmically 
divergent for the hybrid action, can be analyzed in the same way and the 
following bounds are found, 

|r (2) | < C\q\ 2 Log(\q\b); |r (3) | < C\q\Log(\q\b). (98) 

However, as we saw these bounds are not good enough to prove / conver- 
gence due to the discontinuities of the interpolation, and the stricter bounds 
( f44|) are needed. In the continuum, we know that gauge invariance implies 
that the two-point function is of the form, 

r$(q) = (q^-g^q 2 mq) (99) 
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where n(g) is a scalar function that behaves logarithmically for large exter- 
nal momenta. It is clear then that the component = v satisfies the bound 
fl44]). For finite / however, a form like this is not easy to prove for the 
parity-even parts, because we do not have Lorentz invariance. Fortunately, 
the gauge invariant relations of eqs. ( |43] ) are enough. Let us start with the 
two-point function. Using (H) repeatedly we find, 

r (2)/_N r (2) r (2) 

r$M = = E W-§- - E (ioo) 

where the usual definition, 



/ 



(101) 



has been used. Now, the right-hand-side can easily be bounded and we find 
the desired result, 

|rg(<?)| < E l9«lk/j||rg|;i < E \Qa\\qp\Log(\q\), a,/3^». (102) 



The last inequality has made use of eq. (|9q). The P-even three-point func- 
tion can be analysed in the same way. 

Finally, the parity-odd three point function has to be considered. In this 
case we do not have the help of the gauge invariance constraints ( |43| ) , but on 
the other hand, the diagram contains an epsilon tensor. A direct inspection 
of the three-point function reveals the following structure, 

ltl(<Z,</) ~ f ahc e^ 1 {is(q' )s(q 1 )F l (q,q > ) - 2s(q 1 )s{q'p)F 2 {q,q') 

+2s( qi )s(q'p)F 3 (q, q') - s{q^s{q' p )F A (q, q% (103) 

where 

sin 2 (qf/2) sin(q f) 
P S[q) ~ f 



_ ottt ^ s{q) _ ™^ (1Q4) 



and the functions F are lattice integrals with smaller degree of divergence 
than the original integral, for which the following bounds hold, 

l*l(?,?0l < T^Log(\q\), \F 2t3 (q,q')\ < J-Log(\q\) 
\Q\ \l\ 

\F 4 (q,q')\<±-Log(\q\). (105) 

m 
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Using \s(q)\ < C\q\ 2 , \s(q)\ < C\q\, it is then straightforward to check that 

\r i 3»{q,'j')\<\qa\Log(\q\), a + fx, (106) 

This concludes the proof of the bound (H3j) for the P-odd three-point func- 
tion. 
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